Spin dynamics in the frustrated three-dimensional hyperkagome compound GdaGasO^ 
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We present the first neutron inelastic scattering results on the low temperature magnetic state 
of the three dimensional hyperkagome compound Gd3GasOi2 (GGG). GGG is often classified as a 
strongly frustrated system with a manifold of continuously connected states. However, in contrast to 
the expectation of a continuum of gap-less excitations above a spin liquid-like ground state our results 
reveal three distinct inelastic modes found at 0.04(1), 0.12(2) an 0.58(3) meV at 0.06 K. The inelastic 
modes can be attributed to the magnetic ground state with the lowest and highest energy excitations 
showing spatial dependencies indicative of dimerized short range antiferromagnetic correlations. 
Short range correlations, reminiscent of spin liquid-like order, are static within the instrumental 
resolution (50 AteV) and represent 82 % of the spectral weight. Longer range correlations, first 
observed by Petrenko et a/.[T], develop below 0.14 K within the elastic cross section. The short range 
static correlations and dynamic components survive to high temperatures, comparable to the nearest 
neighbor exchange interactions. Our results suggest that the ground state of a three dimensional 
hyperkagome compound differs distinctly from its frustrated counterparts on a pyrochlore lattice 
and reveals a juxtaposition of spin liquid order and strong dimerised coupling. 

PACS numbers: 



In recent years it has become evident that magnetic 
frustration provides an excellent path to novel and ex- 
otic magnetic order [3]? -[5]. Evocative names such as 
spin liquids, spin glasses and spin ice are associated with 
the frustration of magnetic spins. In spin liquids the 
energy scale of interactions between the spins does not 
influence the ordering temperature and due to a mani- 
fold of degenerate states the spins remain fluctuating at 
temperatures much lower than the interaction energies. 
An illustrative example of a spin liquid with a large spin 
value S >> |, a cooperative paramagnet, is Tb2Ti20y. 
Tb2Ti2C>7 remains disordered down to the lowest temper- 
atures and displays a spin liquid state [6] . The excitation 
spectrum reveals, in addition to crystal field excitations 
[7] , a continuum of fluctuating spins that slow down with 
decreasing temperature but remain fluctuating down to 
0.05 K [8 ]. A second example of a cooperative paramag- 
net is the kagome antiferromagnet deuteronium jarosite 
which shows gapless magnetic excitations extending out 
to at least 20 meV with a linear temperature dependence 
of the spin fluctuation rate [9] . These two examples high- 
light the continuum of liquid-like quasielastic scattering 
typically observed in a cooperative paramagnet and con- 
form to the prediction by Moessner et al. that a system 
of classical spins on some frustrated lattices will observe 
a linear temperature dependence of the spin fluctuation 
rate [HEE]. 

Contrary to expectation this work presents, to our 
knowledge, the first inelastic neutron scattering study 
on a three dimensional (3D) hyperkagome structure, 
GdaGa50i2 (GGG) in which we show spin liquid order 
that is concomitant with distinct gapped modes point- 
ing towards singlet-triplet excitations arising from short 
range antiferromagnetic (AF) correlations. These results 



thus shed new light on the spin dynamics of frustrated 
hyperkagome structures. 




FIG. 1: The garnet structure of GGG with two interpene- 
trating hyperkagome lattices showing the interatomic, inter- 
triangular and inter hyperkagome exchange interactions, Ji, 
J2 and J3, respectively. For clarity only the Gd 3+ are shown. 



Structurally, rare earth garnets such as GGG are one 
of very few realizations of a hyperkagome structure, a 3D 
lattice of corner sharing triangles. Another recently cre- 
ated hyperkagome compound is Na^^Os, [12] ■ In fact, 
GGG consists of two interpenetrating hyperkagome lat- 
tices with the triangulated nearest neighbors linked via 
an exchange interaction, Ji, while individual triangles 
are coupled through 2nd nearest neighbor interactions, 
J 2 and the two hyperkagome lattices linked via a third 
nearest neighbor term J3 [T2]. In GGG the magnetic 
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Gd 3+ spins are isotropic (S = 7/2) and are often con- 
sidered as Heisenberg spins due to single ion anisotropy 
of less than 0.04 K [T3]. However the non-negligible 
dipole exchange, D — 0.7 K, could lead to anisotropy 
[T3] , Magnetically GGG shows a Curie- Weiss tempera- 
ture of -2.3 K indicative of AF correlations but does not 
order down to 0.025 K [THUS]. Although indications of 
short range order was hinted at by bulk measurements 
[TBI fTfj] , proof was obtained by neutron diffraction that 
revealed a spin liquid-like ground state down to 0.14 K 
with the development of sharper but not resolution lim- 
ited magnetic diffraction peaks below 0.14 K [JJ in addi- 
tion to the spin-liquid like scattering. Interestingly long 
range magnetic order is achieved via the application of 

only 1 Tesla IIUQII. 

The dynamic nature of GGG has previously been stud- 
ied first via the indirect measurement of /xSR. Two /iSR 
studies confirmed the absence of long range order down 
to 0.025 K, however these studies disagree on the nature 
of the slowing down of the spin fluctuations. In a study 
by Dunsiger et al. [Jj|] , a linear decrease of Gd spin fluc- 
tuations was observed below 1 K which extrapolated to 
8.2 [leV at K while Marshall et al. [20], also determined 
the slowing down of fluctuations but observed a temper- 
ature independent relaxation below 0.2 K. A more di- 
rect study, by Mossbauer spectroscopy , observed the 
fluctuating Gd spins down to 0.027 K with a decrease in 
spin fluctuating rate from 11.9 /zeV at 0.4 K to 0.12 /zeV 
at 0.09 K. Most recently, Ghosh et al. [32] pointed to a 
new dynamical phenomena, in the low temperature phase 
below 0.14 K, in which fluctuating uncompensated mo- 
ments coexist with unsaturated AF order and quantum 
protectorates of defect centered clusters. 

Theoretically Yavorks'kii et al. [53] were able to re- 
produce the spatial correlations of the low temperature 
(T< 0.14 K) phase by taking into account the nearest 
neighbor and the nearly infinite dipole exchange inter- 
actions, J = 0.107 K 13 and D = 0.7 K [24; . Ya- 
vorks"kii et al. showed that J and D are perturbed by 
much smaller exchange interactions Ji and J3 [23j [24] 
and these smaller components dictate the incommensu- 
rate ordering wavevector of the low temperature phase 
T< 0.14 K. The work in this Letter presents inelastic neu- 
tron measurements on polycrystalline GGG from which 
we obtain both spatial and temporal information |25j . 
Neutron time-of-flight measurements were performed at 
the spectrometer IN5 of the Institut Laue-Langevin [5S] . 
IN5 was set up to an incident energy of Ei — 1.94 and 3.27 
meV with average elastic linewidths of 50 and 80 /J,eV, 
respectively, full width at half maximum (FWHM). The 
resolution was determined using a standard incoherent 
scatterer. The temperature dependence of the scattering 
function S(Q,uj) was measured between 0.06 and 9 K. 
The instrumental background was measured using an 
identical empty cell at 2 K and subtracted from the raw 
data. The sample used in this work is that used in the 
previous work of Petrenko et al. [JJ containing 99.98 % 
of the non-absorbing isotope 160 Gd. High resolution neu- 



tron diffraction using D1A of the Institut Laue Langevin, 
A = 1.9 A, was used to determine the upper level of a 
possible disorder on the Ga/Gd sites. The refinement re- 
vealed a fully stochiometric sample, the error of the site 
occupations indicated the upper limit of disorder to less 
than 2%. 

S( Q, w) at 0.06 K is shown in Fig|5] Clearly there 
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FIG. 2: Powder averaged scattering function S(Q, w) of GGG 
at 0.06 K with incident Ei = 1.94 meV. Insert is a cut at 
wavevector transfer = 0.5 A showing the elastic lineshape 
and two low lying excitations with the corresponding fits as 
described in the text. 



is a first inelastic contribution at E = 0.58(3) meV. The 
lower energy contributions are more easily understood by 
taking a cut at constant wavevector transfer, see inset of 
Fig[2] and are well described by a Gaussian elastic line, 
the FWHM of which is fixed by the vanadium standard, 
and two further inelastic contributions at 0.04(1) and 
0.14(2) meV. To avoid confusion the three excited states 
are henceforth named INS1 (0.04 meV), INS2 (0.14 meV) 
and INS3 (0.58 meV). The inelastic scattering cross sec- 
tions can be characterized by the dispersion relation, life- 
time (r ~ r _1 ) and intensity. These parameters can be 
obtained by linking the neutron inelastic magnetic exci- 
tation to the dynamic susceptibility via 

S(Q,u>) = -{n(u) + l)}F 2 (Q) x "(Q,c), (1) 

where F(Q) is a dimcnsionless structure factor that fol- 
lows the magnetic form factor and {n(w)-fl)} is the ther- 
mal population factor. The dynamical susceptibility can 
be further described by a Lorentzian form corresponding 



3 



to an exponential decay of excitations in time, written in 
terms of a damped harmonic oscillator: 



x (Q,w) 
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where O 2 = 



T 2 

9' 



T g is a q-dependent linewidth, cor- 
responding to the FWHM of the peak. Furthermore the 
fit function is convolved with the instrumental resolution 



c 

-ri 

<- 



c 

3 



c 
a) 



0.035 
0.03 

0.025 
0.02 

0.015 
0.01 

0.005 

3 
2 
1 

-1 
-2 



0.06 K-y\ 
0.25 




i 



* 



p. < 

I •« 



0.06 - 0.25 K 



0.6 



1.4 



1.8 



Wavevector transfer (A 1 ) 

FIG. 3: (a) Elastic scattering cross section at 0.06 and 
0.025 K. (b) Difference in elastic scattering cross section be- 
tween 0.06 and 0.25 K. The dashed line corresponds to longish 
ranged ordered peaks observed by Petrenko et al. [T] 



The data has been characterized at each position of 
wavevector transfer. Fig. [3ja) shows the elastic scatter- 
ing observed at 0.06 and 0.25 K which represents 82 % 
of the total scattering. This scattering is reminiscent of 
a spin liquid-like structure factor. In addition, incom- 
mensurate Bragg peaks corresponding to longer range 
correlations develop below 0.14 K within the elastic line. 
The resolution of the elastic line gives an upper limit of 
50 /xeV for the spin fluctuation rate probed. Previous 
/xSR [HI [20] and Mossbauer measurements [? ] indi- 
cate fluctuation rates around 8.2 /ieV which this data 
will not be sensitive to. However, the extra scattering 
that develops below 0.25 K , Fig{3] (b), shows features 
of interest, marked by the dashed lines, representing the 
scattering from longer range static correlations first ob- 
served by Petrenko et al. [Tj and theoretically reproduced 



by Yavorks'kii et al. [53]. Although longer range corre- 
lations exist there is no sign of associated spin waves. It 
is possible that these are too weak to be observed since 
the correlations remain finite on the scale of 100 A pQ. 
The three inelastic peaks, INS1, INS2 and INS3 are all 
dispersionless within the resolution probed. These ex- 
citations do not originate from either local vibrational 
excitations nor crystal field excitations since their de- 
pendence on wavevector transfer neither increases with 
\Q\, as would be expected from local vibrational exci- 
tations, nor follow the Gd 3+ form factor expected for 
crystal field excitations [231 |2"T] . The possibility that the 
higher level excitation is a crystal field excitation affected 
by an internal molecular field can be excluded as this is 
not compatible with the specific heat data [12] [55] . 

FigQa) shows the wavevector dependence of the in- 
tegrated intensity of the three inelastic peaks. The in- 
tensity of INS3 at 0.06 K, integrated in energy across the 
region of interest, reveals spatial correlations correspond- 
ing to the scattering cross section expected from a triplet 
excitation above a ground state of singlet dimers [29] . In 
this case, a dimer can be understood as short range order 
of AF coupled spins within a cluster effectively shielded 
from its neighboring cluster. Magnetic interactions be- 
tween clusters can therefore be neglected. The neutron 
scattering cross section for such a ground state is given 
by 
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where d is the separation between spins, F(Q) is the 
Gd 3+ magnetic form factor [57] and A(T) is a tempera- 
ture scaling factor linked to the canonical partition func- 
tion proportional to the thermal distribution of the sin- 
glet ground state and the triplet excited state, A(T) = 
1/(1 + 3exp(— JodNN /ksT)), with T — temperature, fee 
the Boltzmann constant and JgcINN is the nearest neigh- 
bor exchange energy Jg<inn = JnnS(S + 1) = 1.68-fiT 
[29] . The dashed lines in FigQa) correspond to spatial 
correlations with near neighbor exchange interaction (- 
-), d = 3.7915 A and next nearest exchange interaction 
(-•-) d = 5.7916 A. Clearly the data are well described by 
a model including only near neighbor exchange interac- 
tions. The energy FWHM of this excited state at 60 mK 
is equivalent to the instrumental resolution, thus indi- 
cating strong coupling between dimerized Gd moments. 
The inset of Fig. |4|a) shows the integrated intensity of 
the wavevector transfer of INS1, INS2, obtained with in- 
cident neutron energy Ei = 1.94 meV to optimize resolu- 
tion. The variation of width and position of the peaks do 
not exceed the resolution of the instrument and are thus 
considered as constants. The wavevector dependence of 
the integrated intensity of INS1 follows closely the short 
range order behavior of INS3, the dashed line represents 
the lineshape of Eqn. ??. At low and high wavevec- 
tor transfer the model of short range correlations fails 
indicating that extra terms remain important for a full 
description of the ground state. The integrated intensity 
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of INS2 shows a minima close to the position in reciprocal 
space that corresponds to nearest neighbor correlations 
indicating that the origin of INS2 is very different to that 
of INS1 and INS3. 

The temperature dependence of the normalized inte- 
grated intensities of INS1, INS2 and INS3 is shown in 
Fig|4](b) (top) . The dashed line represents the thermal 
behavior expected from a singlet-triplet excitation, A(T), 
with an exchange interaction JgcINN = 1-68 K, the near 
neighbor exchange interaction obtained by Schiffer et al. 
[IB] , The dashed line follows closely the integrated in- 
tensity of INS3 thus further validating the notion of a 
dimerized short range AF ordered state. Neither INS1 
nor INS2 follow the temperature dependence of INS3. 
The integrated intensity of the INS1 excitation follows 
a trend similar to INS3, albeit with a reduced exchange 
interaction J = 1.3 K, up to 0.6 K but at higher tempera- 
tures does not follow this trend. The integrated intensity 
of INS2 has a maxima at 0.6 K. 

Further information concerning the INS3 excitation is 
revealed in Fig. gb) (bottom) . Unlike INS1 and INS2, 
INS3 shows a strong temperature dependence in its en- 
ergy position and the excitation lifetime. The peak en- 
ergy position can be followed by a power law function 
with parameters j3 — 0.12(1) and falls to zero at T = 
1.67(2) K. The excitation lifetime can be described by 
the algebraic form r = A T 5 with £ = 0.84 ± 0.21 and 
A diverges to the inverse of the instrument resolution at 
0.06 K, see the inset of Fig. Qb) (bottom) . The relevance 
of these parameters becomes clear when reviewing recent 
theoretical work by Robert et al. [30] showing that, in 
contrast to the pyrochlore lattice jTUl [TTJ |3T] , sufficient 
temporal and spatial stiffness in a classical kagome anti- 
ferromagnet can give rise to magnetic excitations corre- 
sponding to acoustic and optical modes in addition to a 
soft mode. 

These excitations depend strongly on the tempera- 
ture regime. At high temperatures (T/J > 0.2 ) only 
a quasielastic signal centered at meV is expected. On 
decreasing T/J from 0.2 to 10 -2 the quasielastic signal 
splits into two excitations, an acoustic mode and a non- 
dispersive soft mode that softens with decreasing temper- 
ature below T/J < 10~ 2 . An optical mode develops at w 
— 2 J for T/J < 5e— 3. Furthermore, the theoretical exci- 
tations are characterized by the temperature dependence 
of their lifetime with an algebraic dependence r = AT^ 
with £ = 0.995 for quasielastic scattering and £ = 0.18 
for inelastic scattering in the regime of cooperative para- 
magnetism, T/J < 0.1. The temperature regime probed 
in this work extends from the partially ordered phase at 
0.06 K into the paramagnetic regime at 1.2 K correspond- 
ing to 0.04<T/J<0.7 (with J = 1.68 K). 

The excitations observed cannot be assigned to acous- 
tic excitations which would remain dispersive and orig- 
inate from Bragg peaks, even with powder averaging. 
However it is possible to assign the non-dispersive gapped 
excitations to the optical or soft modes with a high energy 
mode at ui — 4 J and not 2 J as predicted [3U]. Robert cal- 



culated the excitations for a classical Heisenberg kagome 
AF. It is well known that substantial long range dipole 
exchange interactions play an important role, such a 
mode would therefore be lifted upwards due to further 
exchange interactions and could therefore be observed at 
u) = 4J. Analogous to phonons, a magnetic optical mode 
can arise from a localized perturbation of interactions as 
found in the short range dimerized interactions displayed 
by INS1 and INS3. 
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FIG. 4: (a) Integrated intensity of INS3 as a function of 
wavevector transfer. The dashed lines correspond to a model 
of short range AF correlation with only near neighbor (-) 
or next nearest neighbor correlations (--). The inset figure 
shows the integrated intensities of the two low lying exci- 
tations INS1 and INS2. (b)(top)Temperature dependence 
of the integrated intensity for INS1, INS2 and INS3. The 
dashed line is the temperature dependence expected for a 
singlet to triplet excitation with an exchange interaction 
Jg<inn = JnnS(S + 1) — 1.68K and follows closely the inte- 
grated intensity of INS3. (b) (bottom) Peak position in energy 
of INS3 (dashed line is a fit to the data with a power law corre- 
sponding to a transition temperature, T = 1.67(2) K.), inset: 
T/J dependence of excitation lifetime of INS3. Dashed line is 
the algebraic function r = AT^ with an exponent £ = 0.84 ± 
0.21 



The data presented in this work sheds light on the un- 
usual magnetic ground state of the hyperkagome struc- 
ture GGG. The data indicate that the longer range or- 
der, observed below 0.14 K in previous diffraction work 
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is static and does not impact on the behavior of the 
higher energy spectral density which can be linked to 
the partially ordered state of GGG. 82 % of the scatter- 
ing is static with a fluctuation rate of less than 0.05 /ueV 
incorporating a large component of the spin liquid-like 
structure factor. The remaining spectral weight lies in 
three gapped magnetic excitations two of which can be 



modeled with the spatial dependence of short range AF 
dimer-like correlations. This is highly unusual in a com- 
pound with a ground state manifold in which a contin- 
uum of excitations is expected to be characteristic of the 
dynamic nature of the magnetic ground state, thus lead- 
ing to a new class of magnetic dynamic order for hyper- 
kagome compounds. 
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